
use "A2_Curtailment.dta", clear


**** Appendix Figure A2: Simulated Water Right Curtailment
********************************************************************************
********************************************************************************

colorpalette #fafa6e #2a4858, ipolate(30, HCL) nograph
twoway (rbar y1 y2 short if ID==1 & short!=0, fcolor("`r(p30)'") blwidth(none) barwidth(.09)) ///
(rbar y1 y2 short if ID==2 & short!=0, fcolor("`r(p15)'") blwidth(none) barwidth(.09)) ///
(rbar y1 y2 short if ID==3 & short!=0, fcolor("`r(p20)'") blwidth(none) barwidth(.09)) ///
(rbar y1 y2 short if ID==4 & short!=0, fcolor("`r(p25)'") blwidth(none) barwidth(.09)) ///
(rbar y1 y2 short if ID==5 & short!=0, fcolor("`r(p10)'") blwidth(none) barwidth(.09)) ///
(rbar y1 y2 short if ID==6 & short!=0, fcolor("`r(p4)'") blwidth(none) barwidth(.09)) ///
(rbar y1_cut y2_cut short if ID==1 & short!=0, fcolor(white%70) blwidth(none) barwidth(.09)  yaxis(2)) ///
, ylab(0 "100%" 4845 "90%" 9690 "80%" 14535 "70%" 19380 "60%" 24225 "50%" 29070 "40%" 33915 "30%" 38760 "20%" 43605 "10%" 48450 "0", nogrid) legend(pos(6) col(6) label(1 "Bloomington") label(2 "Gunlock") label(3 "New Santa Clara") label(4 "St. George Clara")  label(5 "Ivins")  label(6 "Lower Gunlock") order(1 4 3 2 6 5)) ytitle("Percent Curtailed") xsize(20) ysize(11) xtitle("Reductions to Total Supply") xlabel(.1 "10%" .2 "20%" .3 "30%" .4 "40%" .5 "50%" .6 "60%" .7 "70%" .8 "80%" .9 "90%", ) ylabel(0 "Senior" 48450 "Junior", axis(2) labsize(med) notick) ytitle("", axis(2))
